The Haldane gap for the 5 = 2 antiferromagnetic Heisenberg chain revisited 



Xiaoqun Wang 1,2 , Shaojing Qin 3 and Lu Yu 4,3 
1 Institut Romand de Recherche Numerique en Physique Des Materiaux (IRRMA), PPH-333, EPFL, CH-1015 Lausanne, 

Switzerland 

2 Max-Planck-Institut fur Physik Koraplexer Systeme, Nothnitzer Str. 38, D-01187 Dresden 
3 Institute of Theoretical Physics, P. O. Box 2735, Beipnq 100080, P.R. China 
4 International Center for Theoretical Physics, P. O. Box 586, 34100 Trieste, Italy 

(February 1, 2008) 

Using the density matrix renormalization group (DMRG) technique, we carry out a large scale 
numerical calculation for the S = 2 antiferromagnetic Heisenberg chain. Performing systematic 
scaling analysis for both the chain length L and the number of optimal states kept in the iterations 
m, the Haldane gap A(2) is estimated accurately as (0.0876 ± 0.0013) J. Our systematic analysis 
for the S — 2 chains not only ends the controversies arising from various DMRG calculations and 
Monte Carlo simulations, but also sheds light on how to obtain reliable results from the DMRG 
calculations for other complicated systems. 



PACS: 75.10.Jm, 75.40Mg 

Since HaldaneEl conjectured that integer-spin 5 an- 
tiferromagnetic Heisenberg (AFH) chains Jaave a gap 
A(5) in the excitation spectrum, much workBll has been 
done on spin one Heisenberg AFH chains. In particu- 
lar, the valence bond solid (VBS) model has been pro- 
posed to provide physical interpretation for this impor- 
tant propertyfl In recent years, very accurate estimates 
of the gap [A(l)n= 0.41049(2) J] have been obtained by 
White and HusejS and by Golinelli et al^i using the_den- 
sity matrix renormalization group (DMRG) methodO and 
exact diagonalization with proper extrapolation, respec- 
tively. Experimental evidence for the gap has also been 
clearly shown on some quasi-one-dimensional S = 1 AFJI 
materials, p±db, as Ni(C 2 H 8 N 2 )pM) 2 (C104) (NENP),EHl 
Y 2 BaNi0 5 py and AgVP 2 S 6 I 3 ] Moreover, nontrivial 
degenerate ground states show up when the chain be- 
comes open. This property is related, to the impurity ef- 
fects observed in Y 2 _ 2; Ca 2 ;BaNi05Jlil As pointed out by 
one of us, the boundary effects and associated impurity 
effects distinguish the low-energy properties of the Hal- 
dane system i;om other quasi-one-dimensional gapped 
spin systems. t3 

The 5 = 2 AFH chain, beyond its own interest, 
is essential for further verification of the Haldane con- 
jecture. In the limit of large 5 the spin chains be- 
come quasi-classical, which implies that the difference 
between the half-integer and integer spin chains, e.g. 
the absence/presence of the Haldane gap, must dimin- 
ish. The 5 — 2 Heisenberg AFH chain is thus the 
first case to check this behavior. Recently there appears 
to be also experimental evidence for a Haldane gap in 
an 5 — 2 system.Ej The material under study is (2,2'- 
bipyridine)trichloromanganese(III), MnCi3(bipy). The 
manganese ions form effective 5 = 2 spins and are cou- 
pled to a quasilincar chain by chlorine ions. The Haldane 
gap has been measured by determining the low temper- 
ature behavior, and it was estimated as 0.07 J. 

For the 5 = 2 antiferromagnetic Heisenberg chain de- 
scribed by 



L 

H = jJ2Si-S i+1 , (1) 

8=1 

numerical calculations are much more elaborate than for 
the 5 = 1 case. There are two reasons. The first one is 
that the value of the Haldane gap for 5 = 2 is roughly 
one order of magnitude smaller than that of 5 = 1 chain. 
Correspondingly, the correlation length in the ground 
state also increases roughly by one order of magnitude. 
Consequently much longer chains are required to reach 
the convergent regime with respect to the chain length. 
In fact, to exclude the finite size effects one has to go up 
to chain lengths of about one thousand for open bound- 
ary conditions. The second reason is that the number of 
degrees of freedom per spin is five for 5 = 2 instead of 
three for 5 = 1, so the exact diagonalization can be car- 
ried only for rather short chains (of the order 10, much 
shorter than the correlation length) . 

Nevertheless, the 5 = 2 Heisenberg AFH chain has 
been numerically studied by quite a few groups using-di£< 
ferent Monte Carlo techniques and DMRG methods ,tZrE2l 
and we have listed results of these calculations chronolog- 
ically in Table I. From these values, on the other hand, 
it is not easy to extract a reliable estimate of the Hal- 
dane gap for the 5 = 2 chain. Recently, Kashurnikov et 
al£3 based on their "Worm" Monte Carlo (WMC) results 
(which are excellent for 5 = 1/2, 1 cases) expressed some 
doubt on the value of A(2) obtained by DMRG calcu- 
lations. This has indeed independently raised a serious 
question concerning both the capability of the DMRG 
method itself and the true value of the Haldane gap for 
the 5 = 2 case. It seems true that analysis based on 
a single DMRG calculation with fixed number of kept 
states m does not lead to a conclusive answer. In this 
paper, we carry out a large scale DMRG calculation 
with up to one thousand and two hundred sites as well 
as a systematic scaling analysis for the number of the 
kept optimal states, up to 400. We will show then that 
A (2) = 0.0876 J± 0.0013 J. In addition, this analysis also 
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sheds some light on how to judge in which case a DMRG 
calculation and the corresponding scaling analysis can be 
trusted. 

TABLE I. Values of the Haldane gap A for S = 2 Heisen- 
berg AFH chain as calculated by various groups using dif- 
ferent techniques. WLMC: World line Monte Carlo method; 
PMC: Projection Monte Carlo Method; LCMC: Loop clus- 
ter Monte Carlo method; WMC: Worm Monte Carlo method. 
BCs: boundary conditions; PB (OB): periodic (open) bound- 
ary conditions, m: the number of states kept in the DMRG 
iterations. SA indicates whether or not a scaling analysis has 
been used to determine the gap value. L indicates the largest 
size of the chains in each case. The numbers in parentheses 
are errors, for the gap at last digit(s), given in those refer- 
ences. In the following table, we set J = 1. 



metho 



idfi-i 



BCs 



SA 



WLM 
DMRGK 
DMRS- 
PMCE-, 
DMRG- 
WLMC- 
DMRG^ 
WLMC- 
DMRG£ 
DMRi 
LCM 1 
WM' 
DMR 



0.08 
0.02 

0.055(15) 
0.05 

0.055(15) 

0.049 

0.085(2) 

0.074(16) 

0.082(3) 

0.085(1) 

0.09(1) 

0.1032(7) 

0.0907(2) 



64 

100 

70 

32 

40 

128 

350 

512 

150 

80 

400 

500 

600 



PB 
OB 
OB 
PB 
OB 
PB 
OB 
PB 
OB 
PB 
PB 
PB 
OB 



No. 
Yes. 
Yes. 
No. 
Yes. 
Yes. 
Yes. 
Yes. 
Yes. 
Yes. 
No. 
No. 
Yes. 



120 
110 

110 

300 

250 
1700 

400 



Let us first recall some related physical .properties. 
Firstxrf all, as emphasized by Schollwock et ac3 and Qin 
et a/J13 for the periodic boundary conditions, the ground 
state is a singlet, whereas for open boundary conditions 
the topological edge excitations make S tot — 0,1,.... 
states non-degenerate at finite lengths. Of course, in 
the thermodynamic limit these states become degenerate, 
and for relatively short chains, the interaction between 
these topological excitations has to .-be carefully taken 
into account in the scaling analysis.EJ Another way of 
removing the "surface" effects is to introduce .different 
spin values and coupling constants at the edgesEj In-our 
numerical calculations following Schollwock et al,cB we 
attached two S = 1 spins to screen the edge spins at 
both ends of an open 5 = 2 chain. This trick is helpful 
in practice for an accurate evaluation of the Haldane gap, 
although the spin gap as a bulk quantity should be inde- 
pendent of the boundary conditions. The main reason is 
that the surface energies with different S 1 ' ' are affected 
by truncations differently than the bulk properties. 

As follows from the one dimensional field theory, the 
leading finite size correction to the Haldane gap with 
open boundary conditions is proportional to the inverse 
of the length squaredfl 



A(m = oo,L) = A(S) 



2A(S)L 2 



where v is the spin wave velocity, and A(5 I ) the value 
of the Haldane gap in the thermodynamic limit. This 
formula shows that the true gap should be a minimum 
of the parabola as a function of the inverse chain length 
i" 1 , if all states are kept (to = oo) and higher order 
terms are negligible. In the DMRG calculations, however, 
the truncation of the Hilbert space leads to deviations 
from the above asymptotic behavior. The deviation is 
not so serious for the spin 1/2 or spin 1 systems, which are 
usually not explicitly emphasized. For more complicated 
systems, like the S — 2 case, the deviation is substantial, 
especially when to is not sufficiently large. 
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FIG. 1. The gap A(m,L) as a function of the number of 
states m kept in the DMRG calculations and the chain length 
L. Each cross for a given m indicates the position of the 
minimum of A(m, L). The solid square denotes the upper 
estimated value Au, while the solid circle denotes the lower 
estimated value Al. The star on the vertical axis indicates 
the recent WMC result .E3 

In Fig. 1, we show the gap A(m, L), as measured 
by the difference between the lowest energy states with 
gtot _ y anc [ gtot _ respectively, as a function of the 
chain length L and the number of states m kept in the 
DMRG iterations. One can immediately see from Fig. 
1 that there exists a minimum as denoted by the cross 
for any given number of kept states to. Therefore, in 
principle one cannot use the extrapolation with respect 
to 1/L in order to get the gap value for the thermody- 
namic limit when m is not sufficiently large. As pointed 
out by Schollwock et al, when the chain length is not 
long enough, the scaling is essentially dominated by the 
1/pLterm . This issue has also been addressed by Qin et 
alE3 Taking this into account one may still obtain a rea- 
sonable estimate of the gap value for the open chains by 
keeping a sufficieptly large number of states, asjipne by 
Schollwock et a/,B Qin et at5 and Schollwock .E3 Their 
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obtained values are comparable with our present results 
shown in Fig.l. 

However, in order to obtain a reliable value for the 
thermodynamic limit, one has to make some further anal- 
ysis as illustrated in Fig. 1. It is obvious that the mini- 
mum of each curve for A(m, L) shifts towards the vertical 
axis, as m increases. In the limit m — + oo the mini- 
mum should be right on the vertical axis, as follows from 
Eq.(2). In fact, the convergence of minimum location 
with respect to m is quite fast. The minimum value of 
A(m = 400, L) is 0.0917385344 at L = 1010B We have 
fitted these positions for m = 120,150,200,250,300,400 
with respect to 1/L. We note that the ^extrapolation 
is performed for m and L simultaneously.!^ The values 
Al and Ajj correspond to two different extrapolations 
using the polynomial fit and the spline fit, respectively. 
We have also checked that the results change very little, 
by leaving out the minimum for m = 120 or by includ- 
ing more minima in the range 150 < m < 400. These 
changes are small compared with the difference between 
Aj/ and A L . 

Therefore, we can safely conclude that the Haldane gap 
is 0.0876 ± 0.0013 for the S = 2 AFH chain, where the 
uncertainty comes from the simultaneous extrapolation 
only. Our analysis is based on the comparison of the 
one dimensional field theory asymptotic and our actual 
DMRG results. 

It is interesting to see how far are the values of the 
Haldane gap obtained by the DMRG method for S = 1 
and S = 2 AFH chains from those given by the follow- 
ing asymptotic formula for large S in the non-linear a 
modekB 

A(5) = aS(S+l)e- ,r V s ( s + 1 > ) for S > 1 (3) 

If we fix the numerical factor a using the value for S = 1, 
then we obtain A(S = 2) = 0.0470, which is only 55% of 
our value. We therefore conclude that the above equation 
is only a semiquantitative, asymptotic formula. 

Before concluding we would like to make two remarks 
regarding the numerical calculations: 

i) It is well-known that Monte Carlo techniques involve 
both statistical errors and systematic errors. In fact, the 
systematic errors must be efficiently eliminated before 
analyzing the statistical errors. For an accurate evalu- 
ation of quantities such as the ground state energy and 
energy gap, it is usually difficult to ensure to which extent 
the configuration dependence is sufficiently eliminated in 
reducing the systematic errors. It is very nice that the 
Worm Monte Carlo simulations have produced excellent 
results for S — 1/2, 1 AFH chains- reproducing even the 
"relativistic" dispersion relationsES for S = 1 case. The 
same method was used to obtain A(2) = 0.1032(7), which 
is quite different from our result. We would like to men- 
tion that the excellent results in S = 1/2, 1 cases cannot 
guarantee that the systematic errors have been efficiently 
excluded for the S = 2 case, since the configuration de- 
pendency involved in these cases is quite different. This 



relevant issue for the Monte Carlo simulations was also 
addressed in Ref. 

ii) The DMRG calculations involve a systematic error, 
too. It originates from the truncation of the renormalizcd 
Hilbert space. In relatively simple applications, like S — 
i and S = 1 Heisenberg chains, these errors are not 
serious, while for more complex systems, such as S = | 
and S = 2 chains, they become rather crucial. With more 
and more applications of DMRG to complicated systems, 
we note that the analysis based on DMRG calculations 
becomes reliable for the thermodynamic limit 

1) when a given set of numerical data follows very 
well the analytic scaling up to a convergent regime 
of scaling variable (usually 1/L); 

2) when various sets of numerical data are computed 
for a number of m like what we have done for the 
S = 2 chain, and a further scaling for m is being 
carefully carried out. 

For the first case, we quote two more examples. One is-fctfj 
a S = 3/2 periodic chain with m=1200 sates keptSo 
The analytical scaling behavior is available from the con- 
formal filed theory, and a careful scaling analysis has 
shown that S = 3/2 spin belongs to the same universality 
class, as S = 1/2 case with central charge c = 1 and the 
scaling exponcpt r) = 1. Another example is an S = 2 
periodic chaipJcJ and the number of states was kept up 
to m = 1700.L3 Again, it has been shown unambiguously 
that S = 2 chain belongs to the same universality class as 
5 = 1 chain. In both cases, the results for the low-lying 
energy excitations and the correlation functions fit the 
scaling functions very nicely in the convergent regime of 
the scaling variables. We can conclude that these DMRG 
results, like the present one for the S = 2 Haldane gap, 
are reliable. 
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